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Abstract 

Presented are the integral solutions to the quantum Knizhnik-Zamolodchikov equations for the 
correlation functions of both the bulk and boundary XXZ models in the anti-ferromagnetic regime. 
The difference equations can be derived from Smirnov-type master equations for correlation functions 
on the basis of the CTM bootstrap. Our integral solutions with an appropriate choice of the integral 
kernel reproduce the formulae previously obtained by using the bosonization of the vertex operators 
of the quantum affine algebra ^^^(sb). 



1 Introduction 

In this paper we address both the bulk and boundary XXZ models in the anti-ferromagnetic regime 
on the basis of Smirnov's axiomatic treatment of massive integrable models Smirnov found three 
axioms that the form factors of the sine-Gordon model should satisfy The said three axioms consist 
of the S-matrix symmetry, the cyclicity condition, and the annihilation pole condition. Smirnov also 
pointed out in that the first two axioms imply the quantum KZ equations |^ of level for form 
factors. 

We have a similar story for correlation functions under the scenario based on the CTM (corner 
transfer matrix) bootstrap approach Q . As was shown in that correlation functions of bulk massive 
integrable models such as the XYZ model should satisfy the three relations: the R-matrix symmetry, 
the cyclicity condition and the normalization of correlation Junctions. Let us call those three relations 
Smirnov-type master equations. The first two relations imply the quantum Knizhnik-Zamolodchikov 
(KZ) equations |^ of level —4. The boundary analogue of Smirnov-type master equations were also 
obtained in |^ for the boundary XYZ model. In the boundary model case, the second cyclicity condition 
should be replaced by reflection properties. 
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We start from Smirnov-type master equations for correlation functions of the bulk/boundary XXZ 
model in the antiferromagnetic regime. We will make the Ansatz in the subsequent sections that 
the correlation functions are expressed in terms of integral transform of certain rational functions with 
appropriate integral kernels. We determine the said rational functions such that the i?-matrix symmetry 
holds. We also find the transformation properties and recursion relations for the integral kernels by 
imposing the cyclicity conditions (resp. reflection properties) and the normalization conditions for the 
bulk (resp. boundary) XXZ model. 

The above mentioned Ansatz was already made in |0, H, ^ , where we constructed the form factors 
of the anti-ferromagnetic XXZ model by solving Smirnov's three axioms. We will see below that the 
integral formulae for correlation functions and those for form factors have the very similar structure. 
The difference between those formulae mainly results from the choice of rational functions, while the 
integral kernels are essentially the same. For the case of form factors, the freedom of solutions to the 
quantum KZ equation corresponds to that of local operators of the XXZ model. 

In the framework of the representation theory of the affine quantum algebra jlOj the space of states 
of the anti-ferromagnetic XXZ model can be identified with the tensor product of irreducible infinite 
dimensional modules, one is level 1 highest and the other is level —1 lowest modules of C/g(s[2) with 
< —q < 1. Here the restriction < ~q < 1 is required by the anti-ferromagnetic condition. The 
correlation functions of the anti-ferromagnetic XXZ model are given by the trace of the product of the 
type I vertex operators on the space of physical states jl^. Furthermore, the integral formulae for the 
correlation functions of the XXZ model was also obtained in [0 by using the bosonization of vertex 
operators of Uq{5l2). The vertex operator approach were also applied to obtain the integral formulae 
for correlation functions of the anti- ferromagnetic XXZ model with a boundary 

One of advantages of the vertex operator approach is that correlation functions in this framework 
naturally satisfy the quantum KZ equations by the construction. Another advantage is that it is also 
relevant to massless models. In the integral formulae were conjectured by solving the quantum KZ 
equations of level —4 for correlation functions of the massless XXZ model, which has C/g(5[2)-symmetry 
with \q\ = 1. This was based on the CTM bootstrap approach [Q, which is not directly applicable 
to the massless regime. In JTi] those conjectured formulae were reproduced by using the bosonized 
vertex operators. The same bosonization method was used to obtain the corresponding formulae for the 
massless XXZ model with a boundary in 

From the view point of solving the quantum KZ equations, the trace construction allows us to chose 
the level arbitrarily, but gives only a canonical solution. On the other hand, Smirnov-type integral 
formulae based on the form factor bootstrap approach have the freedom of solutions but only solve 
the level quantum KZ equations. As for correlation functions, Smirnov-type integral formulae give 
solutions to the level —4 equations. The relations among those different type integral formulae were 
discussed by Nakayashiki et al Furthermore, Nakayashiki proved that the trace functions of type I 
and type II vertex operators associated with the C/g(s[„)-symmetric model give a basis of the solution 
space of the quantum KZ equation at arbitrary level . The CTM bootstrap approach was developed 
for A^^^i-symmetric elliptic models without and with a boundary, in order to derive the quantum KZ 
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equations of level -2n |T8[|9[ |o). 

The rest of the present paper is organized as follows. In section 2 we formulate the bulk XXZ model 
in the anti-ferromagnetic regime, and present Smirnov-type integral formulae for correlation functions. 
In section 3 we prove that the formulae given in section 2 actually solve the quantum KZ equations. In 
section 4 we further present Smirnov-type integral formulae for correlation functions of the boundary 
XXZ model. In section 5 we give some concluding remarks. In Appendix A we give a simple proof of 
Proposition |[ 



2 Integral formulae for the bulk XXZ model 

2.1 The Hamiltonian and the i?-matrix 

In this section we consider the XXZ spin chain in an infinite lattice 

Hxxz = E(^J+i'^J + + ^"hi'^D- (2-1) 

iez 

Here (t|, ctJ and cr| denote the standard Pauli matrices acting on j-th site, and we restrict ourselves to 
the anti- ferromagnetic regime: A < — 1. The XXZ Hamiltonian Hxxz commutes with Uqish), where 
A ~ (q + q^^)/2 and —1 < q < 0. For lator convenience we also introduce the positive parameter x = —q 
such that < a; < 1. Let V — Cv^ + Cu_ be a vector representation of C/g(5l2). Then the Hamiltonian 



( prl ) formally acts on ^ . . . (g, V V ® ■ ■ ■ . In fig] the space of states V^°° was identified with 

the tensor product of level 1 highest and level —1 lowest representations of {/^(sb). 
Let us introduce the i?-matrix of the six vertex model, where R{C) G End(V^ (E> V): 

R{C)v,,(^v,,= V v,,^v,,RiCTe\i, i?(C) = ^i?(C), (2.2) 



where 



/'4. 4\ f 2 —1 . 4\ 

^(C)-C . 4 1 4. 4^ ' (a;pi,--- ,p„)oo= ll(l-aPi^---p„"), (2.3) 



ki>0 



and z = (^^. The nonzero entries are given by 



R(C)Xz=R{0--t = biO = ^^^^, (2.4) 



i?(C)i; = R{0+t = ciC) 



(i-x2)c 



1 - X'^C 

Assume that the spectral parameter C and the constant x lie in the principal regime: < x < <!• 
As is well known, the XXZ Hamiltonian in the anti-ferromagnetic regime can be obtained from the 
transfer matrix for the six vertex model in the principal regime, by taking logalithmic derivative with 
respect to the spectral parameter 
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The main properties of the i?-matrix are the Yang-Baxter equation 

i?12(Cl/C2)i?13(Cl/C3)i?23(C2/C3) = i?23(C2/C3)i?13(Cl/C3)i?12(Cl/C2), (2.5) 

where the subscript of the _R- matrix denotes the spaces on which R nontrivially acts; the initial condition 

i?(l)=P; (2.6) 

the unitarity relation 

i?i2(Ci/C2)i?2i(C2/Ci) = 1; (2.7) 

the Z2-parity 

i?i2(-C) = -a^i?i2(C)^^; (2.8) 

and the crossing symmetries 

^2\(C2/Cl) = a-,Rl2{x-\l/C2)<yf, i?2\(C2/Cl) - (»a^)i?12(-X-'Cl/C2)(*^D- (2.9) 



The properties (2.6-2.9) hold if the normalization factor of the _R-matrix satisfies the following relations: 

<CHC') = 1, <C)<exO = ebiO, (2.10) 

where e = ±. Under this normalization the partition function per lattice site is equal to unity in the 
thermodynamic limit p^ . 

2.2 Correlation functions and difference equations 

Let us introduce the V"®^"-valued correlation functions 

Gi")(Cl,--- ,C2„)= J2 ^e,^---^Ve,„Gi"\Cl,--- X2nr'-''\ (^ = ±)- (2.11) 

e,-=± 
ei + -- + e2„=0 

Here, we restrict G^,^\C) to the 'total spin-0' subspace of V®'^'^. In the framework of the representation 
theory of C/^(s[2), the correlation function ( ^.11 ) gives the expectation value of the local operator of the 



form 



O ^ e'^^\ ■ ■ ■ E'fl , 



where is the matrix unit on the j-th site, by speciahzing the spectral parameters as follows: 



(0), = Gi")(x-iC,--- ,a;-'C,C,--- ,C)-^"--^^^'^-"". (2.12) 

In what follows we often use the abbreviations: (C) = (Ci, • • ■ ,C2n)7 (C') = (Cij ■ ■ ■ iC2n-i), (C") = 
(Ci, ■ • ■ ,C2«-2), {z) = (zi, ■ • ■Z2n), {z') = (zi , . ■ • Z2„_i ) ; and (e) = (ei---e2n)> {^') = (ei • ■ ■ £2n-i), 
(e") ~ (ei ••■ £211-2). On the basis of the CTM (corner transfer matrix) bootstrap approach, the 
correlation functions satisfy the following three conditions B: 



1. R-matrix symmetry 

P„+iGi")(... ,0+1,0, •••) =i?,,+i(0/0+i)G'i")(--- ,0,0+1, •••) (l^j^2n-l), (2.13) 



where P{x ® y) = yE) x. 
2. Cyclicity 



3. Normalization 



Pl2 ■ • • P2„-12„Gi"HC',a;2C2„) = (TGi")(C2„,C')- 



Gi")(C",C2n-l,C2n)|f,„=..-iC.„-, = Gt''\n ® (e = ±), 

where = D4. (g) u_ + eu_ . 

These three conditions can be componentwisely recast as follows: 

= ^(o/o+imriG^'H--- ,0,0+1, •••)-^^^^+^-- 



(2.14) 



(2.15) 



(2.16) 



Gi"H0,a;'C2„)^'^^" = aG^")(C2«,0)^="^'- 

Combining ( |2l^ ) and (|1|) we obtain another expression of the normahzation condition: 
^ ,(i-)/2G(")(0',C2„-i,e<2„-i)^"^-^ = G-r^HC")^", 

where we also use 

R{ex-^) 



(2.17) 
(2.18) 
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Since the i?-matrix preserves the 'total spin' of the correlation function, we can restrict Gi"^(C) to 
an element of the 'total spin-0' subspace, when wc consider the equations (2.13-0^). Note that the 
first two conditions imply the difference equation of the quantum KZ type H] of level —4: 



X-Rj2n(0/C2n) ' ' ' R-j j + l{C,j / Q + l)Ga'^\() , 

where Tj is the shift operator such that 

r,-F(0 = i^(0,-- - ,x-2o,--- ,C2„), 

for any 2n-variate function F. When the definition of Tj is replaced by 

T,F(0 = F(Ci,--- ,C2„), 



(2.19) 



(2.20) 
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and all the arguments x '^Q/Ck in the first line of the RHS of ( 2.ig| ) are also replaced by x^^'^Q/(k, the 



difference equation (2.19) is called the quantum KZ equations of level /. 



Remark. As a result of Z2-symmetry of i?-matrix, there are two ground states for the anti- 
ferromagnetic XXZ model. Let us specify the two ground states by i = 0, 1, and denote the correlation 
function on the i-th ground state by G,["'' (C) . The CTM bootstrap approach suggest that both Gg""* (C) 
and Gi"^(C) wiU appear in the cyclicity condition as follows: 



Pl2 ■ ■ ■ -P2n-12nG-"''(Cl, ' ■ ' , C2n-l,X^C2n) — G^"''j (^2^ , Cl : " ' ' > C2n-l) 



Thus we introduce Ga^\Q = Gq''{(^) + a'Gi'''{() such that the second equation (2.14) involves only 
Gi")(C). 



2.3 Integral formulae 

Set 



l^j<fesJ2n 

Here c„ is a constant which will be determined below, and the function g{z) has the properties 



The explicit form of g[z) is as follows: 



{x^z;x'^,x^)oo{x'^z ^ix^jx"^) 
{x^z;x^,x^)oo{x^z-^;x^,x'^) 



(2.21) 



(2.22) 



(2.23) 



Thanks to (2.22) the first two equations (2.16- pTT7| ) are rephrased in terms of G^"''(C) and i?(C) as 
follows 



^ ^ 3 ■ J + 1 



(2.24) 



2n-l . 

Gl"^(c',x^C2«)^'-" = Gi"\c2„,c'r"^' n T^- 

j=i ^1 



7=r(")/ 



(2.25) 



In order to present our integral formulae for G^ (C) let us prepare some notation. Let 

A := {a\ea > 0, 1 a ^ 2n}. (2.26) 

Then the number of elements of A is equal to n, because we consider only the 'total-spin-0' case. We 
often use the abbreviations (w) — {wai , • ■ ■ , Wa„), [w') = [wai , ■ ' ■ j Wa,^_-i) and {w") = (wai , ■ ■ ■ , Wa„_2) 
for Gj G A such that ai < • ■ • < a„. Let us define the following ^^^"-valued meromorphic function 



iw\o= J2 ^e,®---®v,,^Q^-\w\cr 



(2.27) 
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where 

f a-l 2n 



q'-'^hmcy = n n(^^j- - n - / n (^^'^^ - ^^-)- (2.28) 

aeA \J = 1 j=a+l J I 

Note that the Qn{o\0Y in (l3), the massless analogue of Q*^"-* (wjC)^, can be obtained by substituting 
Cj = e~^^\ X = e"^" , and Wa = e"^"^"" into (2.28), up to a trivial factor. 
We wish to find integral formulae of the form 

Here, the kernel has the form 

^'^-\w\0=d^-\w\C) n 11^(7^)' (2.30) 



a£A j = l 

where 

^(^) = I . 4w^ -1. 4^ ■ (2-31) 

yXZ^X joo\XZ ,X joo 

For the function ^^^\w\C,) we assume that 

• it is anti-symmetric and holomorphic in the Wa G C\{0}, 

• it is symmetric and meromorphic in the Q G C\{0}, 

• it has the two transformation properties 

a(n) / I*/ 2 a \ 2n — 1 ^ 2n 

v<y {MC C2n) -r-r -pr (2™ _4« TT TT /^-i qo^ 
,(n), = ^11^11^='^^ ' (2.32) 



(uAx«^ ^ ,2„-rT^^^-a„-|-r^, (2.33) 



• it satisfies the following recursion relation 

d^a\w' ,X~^ Z2n-l\C," , C,2n~l, <iX~^ C,2n-l) 



2n-2n 



Here 

mGZ 

and we also fix the constant c„ as follows: 

(_l)"(«+l)/2 (^2.^4^^4)2„ 



1 n n W-^Qx^{xWa/z2n-l)- (2.34) 



a.n(«-l)/2(^2. a:2)^"-l)/2 (3,8. x^) 



00 

2n ■ 

00 



(2.35) 
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The function t?^"-* {w\(^) is otherwise arbitrary, and the choice of ^ji""* {w\() corresponds to that of solutions 
to (2.13- pTT5 ). The transformation property of ii}'^\w\C) iniphes 



aeA 
2n 

n " 



XZ2n 



X^Z2n 



XWa 



2n 

n 



C2 



a . T ^7 



(2.36) 
(2.37) 



The integrand may have poles at 



X 



±(l+4fc) , 



x'^Wf, 
X^'^Wb 



(l<js^2n,fceZ^o), 
{b<a), 
ib>a). 

We choose the integration contour Ca with respect to Wa {a G A) such that Ca is along a simple closed 
curve oriented anti-clockwise, and encircles the points x^'^'^''zj{l ^ j ^ 2n,k £ Z^o) f^nd x'^Wb {b < a) 
but not x^^^'^^ Zj{l ^ j ^ 2n, fc G Z^o) nor x~'^Wb [b > a). Thus the contour Ca actually depends 



on Zj's besides a, so that it precisely should be denoted by C'a{z) — Ca{zi,- 



, Z2r, 



The LHS of 



(2.14) refers to the analytic continuation with respect to C2n- Nevertheless, once we restrict ourselves 



to the principal regime < a; < C "'^ < 1; we can tune all Ca's to be the common integration contour 



C 



because of the inequality xzj < x ^ < x 



3 Main proposition 

In this section we prove the following proposition in the bulk anti-ferromagnetic XXZ model case: 

Proposition 1 Assume the properties of the function •d')j'^ below (2.31) and the integration contour Ca 
below Then the integral formulae ^Jlj with (|!||, ^Idj , \2.3l\j solves the 

three equations ( 2.1^ - 2.15 ). 



3.1 Proof of the i?-matrix symmetry 



Let us first prove (2.13), or equivalently (2.24), in a componentwise way. 



Suppose that j, j + 1 ^ A. Then the relation ( p. 241 ) holds, because the integrand {() is symmetric 
with respect to and Ci+i- 



When j ^ A and j + 1 G A the relation (2.24) reduces to 



1 Cj+iiCji ■ ■■)' 



(3.1) 



- 5(C,/0+i)Gr(--- ,0,0+1, •••)•■-+■■ 
Note that the set of the integration variables in the second term of the RHS is different from the other 
terms. Since Wo's are integration variables, we can replace Wj in the second term of the RHS by Wj^i. 
After that, the relation (|3.lD follows from the equality of the integrands. In this step we use 



Vj+l 



x{zj 



Zj+lj 



Z j -\-\ X Zj 



[xZj 



{l-x^)z, 



Z j -|- 1 X Zj 
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The case j & A and j + 1 ^ A can be proved in a similar way to the previous case. Here we use 

X{Z^ - Zj + i) (1 - X^)Zj + i 

xz.+i - Wo = ^ — (xwo — H ^ — — w.). 

Zj+l-X^Zj 

Finally consider the case j, j + 1 G A. Since is antisymmetric with respect to Wa's, only 

the antisymmetric part of Q^"'\w\C,)"'^^"' with respect to Wj and Wj+i gives non-zero contribution for 



the integral. Thus the relation ( 2.24 ) in this case follows from the fact that the following combination 



{xWj — Zj+i){xZj — Wj + l) {xWj+i — Zj+i){xZj — Wj) 
X~^Wj + l — XWj X~^Wj — XWj^i 

is symmetric with respect to Zj and Zj+i. 
3.2 Proof of the cyclicity 

Let 2n ^ A. When the integral ( |2.29[ ) is analytically continued from (211 to x'^(2n, the points 



move to the points 



9 5 —1 —5 —9 

7 X Z2m X Z2ni XZ2n: X Z2n^ X Z2m X Z2m 



13 9 5 3 —1 —5 

, X Z2n^ X Z2m X Z2n^ X Z2ni X Z2n^ X Z2m 



In the LHS of (2.25) the point x^Z2n and consequently xz2n are outside the integral contour = 
Ca{z' ,x'^Z2n)- Nevertheless, we can deform to the original one Ca = Ca{z) without crossing any 
poles. That is because the factor 

Y[ {xWa - x'^Z2n), 

contained in Q'''^\w\C,' , x'^ C,2nY ~ cancels the singularity at Wa = x'^Z2n- Thus the integral contours for 
both sides of ( ^.25D coincide. 



Furthermore, by using ( 2.36| ) we obtain 

2n 



which implies that the integrands of both sides of ( ^.25 ) coincide and therefore the relation (2.25) holds 
when e2n < 0. 

Next let 2n G A. In this case we also make the rescale of variable W2n ^ x^W2n in the LHS of ( 2.25 ). 



Then the integral contur with respect to Wa (a G A\{2n}) will be = Ca{z' , x'^Z2n), and the other 
one will be C = C2n(x^'^z' , Z2n)- For a £ ^\{2n}, we can deform the contour to the original Ca 
without crossing any poles, as the same reason as before. The integral contour C encircles x^^'^'^'^Zj 
(1 j =^ 2n - 1, fc e Z^o) and x^+^''z2n {k G Z^o), but not a;-5-4fe^^. (i <; j <; 2n - 1, fc G Z^o) nor 
a;"-^~'"^Z2„ (fc G Z^o)- Since Q'^"'^ {w\(' , x'^ (2nY ^ contains the factor 

2n-l _ 

X^C2n TT {xz-i -x'^W2n) TT T"! ~ ' (^-2) 

J--^ LL x~^x'^W2n - XWa 



9 



the pole at W2n — x^^Z2n (1 ^ j ^ 2rj — 1) disappears. Thus we can deform the contour C to the 
original one C2n = C2niz) without crossing any poles. Thus the integral contours in both sides of ( p. 25 ) 
coincide. 



Replace the integral variables such that {w',W2n) * {w2n,w') in the RHS of (|2.25|), and compare 



the integrands of both sides. From (2.36, 2.37) and the antisymmetric property of vE'ct"'' with respect to 
Wo's, we have 

,-T.(n) / I A 1^/ 9> \ 2n— 1 ^ 

^(W ,a: W2„|C ,X Q2n) ^ [^l^-^x^^ TT " '^^^n T-j- C2« Zj - XW2n ,g 

^^a\'W2n,w'\C) oel ^^^2n " a;Wa Q X^W2n - Z^' 



Note that (3.2) describes all W2n- and Z2n-dependence of Q''^^\w\C,' ^ x"^ C,2nY ^ ■ The product of (3.2) and 



^ is equal to 

2n-l ^ 

y. TT S2n , N TT XZ2n ~ '^a 

^2" 11 ^{xW2-.-Z,) W —- 



which implies that the integrands of both sides of ( ^.25 ) coincide and therefore the relation (2.25) holds 



when £2n > 0. 



3.3 Proof of the normalization condition 



Let us prove ( 2.15| ), or equivalently ( 2.18| ). The factor g{z2n~i/ Z2n) hase a zero at (^2n = ^x ^C2n-i- On 



the other hand, the two points xz2n and x~^Z2n-i are required to locate opposite sides of the integral 
contour Ca, so that a pinching may occur as (^2n — * ^x^^C,2n-i- If there is no pinching the correlation 
function G^^\c) vanishes at C,2n — ix^^C,2n-i- 

Suppose 2n — l,2n ^ A. In this case the factor {xwa ~ ^2n-i) contained in (5*-"-'(w|C)^ cancel 
the poles of ip(wa/ Z2n-i) at Wa — x^^Z2n-i, and therefore no pinching occurs as (2n x^^(2n-i- Thus 



the condition ( p.l8[ ) holds if 2n - 1, 2n ^ A. 

When 2n — l,2n G A, no pinching occurs for the integral contour for Ca {a G A\{2n — l,2n}), 
as the same reason as before. Concerning the integral with respect to W2n-i and W2m there is no 
singularities at i«2r!-i = xz2n and W2n = xz2n and consequently no pinching occurs actually. In order 
to see such vanishing singularities, let us first evaluate the residue at W2n = xz2n- Because of the 
antisymmetric property of 1?^"^ with respect to WaS, the zero of i?i"^(ii;", W2n_i, a^22n|C) at u'2n-i = 
xz2n cancels the poles of ip{w2n~i/z2n) at the same point. Furthermore, the poles of '>jj{w2n/ Z2n) 
and ip{'W2n/ Z2n-i) at W2n ~ xz2n and Z2n = x^^Z2n-i are canceled by the zeros of g{z2ii-i/ Z2n) and 
Q^'^\w" ,W2n-i, xz2n\CY at the same points. The latter cancellation can be shown as follows. Note 
that the integral is invariant as we replace Q^'^\w" , W2n-i, xz2n\Cy '^'^ by its antisymmetric part with 
respect to t«2n-i and xz2n- When Z2n — x~'^Z2n-i the antisymmetric part contains the following 
vanishing factor: 

{xW2n-l - Z2n){xZ2 n—1 XZ2n) {xXZ2n ^2n)('^^2n — 1 ^2n — l) 
X-^XZ2n - XW2n-l X-'^W2n-l " XXZ2n 

Thus there is no singularity at W2n = xz2n as C,2n ^ ^x^^C,2n-i- The same thing at W2n-i = xz2n can 
be easily shown in the same way. Hence the condition (ET8) is verified when 2ri — 1, 2?i G A. 
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Next let 2n — 1 ^ A and 2n e A. In this case there is no pinching for the integrals with respect to 
Wa (a G A\{2n}). Let C denote the integral contour with respect to W2n such that C encircles the same 
points as C2n does but xz2n- Note that no pinching occurs with respect to the integral along C because 
both the two points x~^Z2n-i and xz2n lie outside the contour C. Thus the integral with respect to 
W2n along the contour C2n can be replaced by the residue at W2n = xz2n- 

In order to evaluate the residue, the following formulae are useful: 

(^6. „4\ ( 4. 4 4\2 



lim I 



Z2n-1 
Z2n 



XZ2n 



Z2n-iJ (x4;x4)oo (a;2;a;4,a;4) 



„4. ^4 ^4\2 

/CO 

2. ^4 T,4^2~' 

oo 



(3.4) 



Res ^ 



W2n \Z2nJ (x^;^^) 



1 



(3.5) 



g{z)g{x'^z)%l;{xz) 



1 — 



(3.6) 



/ X w \ , , — 
V" - (2: - a;w) = J—- — 

\Z) \ Z J [XW/ Z\X'^)oo\XZ/W\X'^} 



Using these the condition ( 2.18 ) with £2n-i < and e2n > reduces to 



Cn-i^t~'\w'\C') = c„i?^"Hi«',a;-i^2n-i|C",C2„-i,ea;-iC2n-i) 



□ (") 



X {-l^x- 



(^S. „4 ^4^2 2n-2 



.2. ^4 ^4A2 2n- 



(x2;a;4,x4) 



1 n n o 



Wa{x'^;x'^) 



1 = 1 aeA 



ea;2(xWa/z2n-l)' 



(3.7) 



which is valid under the assumption of ( 2.34 ) and ( 2.35 ). 

When 2n — 1 £ A and 2n ^ A, the only difference from the previous case is that the rational function 
Q^"''(w|C)'^ contains the factor 



C2n-l(a;W2n-l " Z2n)\w2n-i=xZ2„ — {^'^ ~ ^)C2n-lZ2n, 

in the present case, while the corresponding factor in the previous case is 



C2n(a;Z2n-l " W2n) 



W2n=XZ2„ 



i:C2n{z2n-l — Z2n)- 



(3.8) 



(3.9) 



Since (3^) is equal to e times (3^) as C,2n — ^x ^C2n-i, the condition ( 2.18 ) with £271-1 > and 62,1 < 



follows from the previous case. After all, the condition ( 2.18| ) was componentwisely proved 



3.4 Non-trivial theta function 



In subsections 3.1-3.3 we proved Proposition The key point in proving (2.16) was the i?-matrix 
symmetry of the rational function Q^"^ (wjC)^. We observed that the other two conditions ( p. 17 ) and 



(2.18) hold under the assumption of the transformation properties and the recursion relation of da 



(") 



But actually, you can easily find ( 2.32 ), ( 2.33 ) and ( 2.34 ) by imposing ( 2.25 ) for £2,1 < 0, ( p.25| ) for 
£2n > and ( ^.17 ), respectively. 
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In this way we obtained the integral solutions to Smirnov-type master equations for correlation 
functions of the anti- ferromagnetic XXZ model, with the freedom of the choice of i^i"'. Now we wish to 
present an example of satisfying all the properties given below (2.31). 



2n 



2n 



i)^-\w\o = e,. -a n n n n '©.^km). 



(3.10) 



a<b 



You can easily see that ( 3.10| ) satisfies all the properties of symmetry with respect to Cj's, antisymmetry 



with respect to WaS, ( ^.32| ), ( ^.33| ) and ( 2.34 ). Note that this example of the function i^^""* (wjC) Jl 
coincides with the theta function for the anti-ferromagnetic XXZ model form factors in the 'total 
spin-0' sector, up to a constant. We therefore conclude that Smirnov-typc integral solutions for both 
form factors and correlation functions hold the essentially common integral kernel, and that the only 
difference among both formulae results from the choice of the rational functions. 

Taking into account that x = in while we set X ~ —q, you can easily confirm that our integral 
formulae with the choice of the function i^^"-* (3.10) reproduce the ones obtained in the trace construction 
up to a trivial constant. Here, the following formulae are useful in order to compare eq. (8.21) in 
[0 with our solutions: 

2n 

aeA j=l 



\2n n 



and 



-(e,2(-u)±e,2(u)) = 



We also notice that there is actually no poles at Wa 
we fix d^a^ to (3.10) because of the factor Q^2(wa/wb)- 



e^8(-a;2u2), 

uQj.a[ — x^V?). 

- x'^wi, (b < a) and Wa = x^^Wf, {b > a) when 



4 Integral formulae for the boundary XXZ model 

4.1 The Hamiltoninan and the X-matrix 

In this section we consider the XXZ spin chain in a half-infinite lattice 

oo 



2 



where 



A = -" + "'\ /z=l^i±I. (4.2) 
2 ' 4a; 1 - r ^ ' 

We again restrict ourselves to the anti- ferromagnetic regime: < x < 1. We also restrict the discussion 
below to the case < < |r| < < 1 such that h > 0. Since the quantity h denotes the magnegic 
field at the boundary site, the Z2-symmetry is broken by the boundary term. In what follows we assume 
that there is only one ground state for the anti-ferromagnetic XXZ model with a boundary. In the 
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presence of the positive magnetic field at the boundary, we should fix the unique ground state to \0)b, 
in the terminology of Ref. . 



As done in the bulk case, the Hamiltonian (4.i) can be obtained from the transfer matrix for the 



six vertex model on a half-infinite lattice, by taking logarithmic derivative with respect to the spectral 
parameter C. The boundary interaction in the boundary six vertex model is specified by the following 
diagonal reflection matrix |2^, 

KiOv, ^ ^-'K{Oi. KiC) = jT^KiC). (4.3) 
Here the scalar function /(z;r) is given by 

and 



{x^rz-x^)^' ^'^^ >^ {x^z^-x^). 



Mz;r) ^ z'2z/r ^ M^') ^ zz^z^r ■ (4-4) 

The non-zero entries are 



1 — 'p'y 

K{0+^k+{z)^-—, X(C)I -fc-(z) = l. (4.5) 
We also introduce another iiT-matrix K'{C): 

K'iOv, = Ve'K{x-\)-_i. 

e' = ± 

The explicit form of K'{C) is given as follows: 

K'iO = 777^^(C), /'(.; r) = x-^z ^i(--V) ^^(^^^-;V) ^ (4 

j'[z;r) ipi{x ^z;r) ip2{x ^z;r) 

and the non-zero entries are 

K^iOt = = ■ = k'_{z) = 1. (4.7) 

K+(a; ^z) 

The main properties of the ii'-matrix are the reflection (boundary Yang-Baxter) equation 

K2{Z2)R21{ZI + Z2)Ki{zi)Ri2{zi " Z2) = i?2l(2l " Z2)i^l (^1 )i?12 (^1 + 22)^^2(^2). (4.8) 

the initial condition 

K{l) = I- (4.9) 



and the unitarity relation 



the Z2-parity 



K{OK{C') - 1; (4.10) 



K{-O^K{0; (4.11) 
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and the boundary crossing symmetry 

X(C)?;= J2 RiC%;illKiex-\)i (e = ±). (4.12) 

Under this normalization the partition function per boundary lattice site is equal to unity in the ther- 
modynamic limit p^ . 

4.2 Correlation functions and difference equations 

Let us introduce the y^^"-valued correlation functions 

Gi"^(Ci,--- ,C2«)= ®•■•®^'e.„Gi")(Cl,•■• ,C2„r^-"^". (4.13) 

£l + ---+e2n=0 

Here, we restrict g[^\c) to the 'total spin-0' subspace of V"®^". The subscript b refers to quantities in 
the boundary model. In the bulk model, we considered two correlation functions G^^\Q for each n. 
On the contrary, because of the Z2-symmetry breakdown, we will consider only one correlation function 
g["^(C) for each n in the boundary model. As in the bulk model case, by specializing the spectral 
parameters the correlation function ( [1.13| ) gives the expectation value of the local operator as follows: 



(^S' •••^ill)fc = Gi")(x-iC,--- ,2:-iC,C,--- ,C)-^"----^^^'^---^". (4.14) 



In this section we often use the abbreviations ((^) = (C2, • • ■ , C2ri), and (e) = {s2, ■ • ■ , £2n)- For fixed 
indices ai, • ■ • , a„ such that ai < ■■ ■ < an, we also use the abbreviation (w) = {wa2, ■ • ■ jf^a„)- On 
the basis of the boundary CTM bootstrap approach, the correlation functions satisfy the following four 
conditions ||]: 
1. R-matrix symmetry 

P„-+iGi")(. • • , 0+1,0, • • • ) = i?„-+i(0/0+i)Gi"'(- • • , 0:0+1, • • • ) (1 ^ J 2n - 1). (4.15) 



2. Reflection properties 



and 



if2n(C2„)Gi")(0,C2n) = G^") (C , ) ; (4-16) 



K[{C,)Gf\c,^\0 = Gf\x-\,X)- (4.17) 

3. Normalization 

Gi"HC",C2„-i,C2„)lc„.=..-c.„-: = Gi"-^HC")®". (e-±). (4.18) 

Since the i?-matrix and if -matrix preserves the 'total spin' of the correlation function, we can regard 
G{;"''(C) as the element of the 'total spin-0' subspace, when we consider the equations ( 4.15 - 4.18 ). Note 
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that the first three equations imply the boundary analogue of the quantum KZ difference equation: 

X i?i,(CiO) • • • Rj-i,{Q-iQRj+io iCj+iQ) ■ ■ ■ i?2n,(C2„0) (4-19) 



where Tj is the shift operator defined by (2.20) 



4.3 Integral formulae 

Set 



2n 



G\r\0 = n CM^j/zk)9{zjZk) n Mzr.r)9biz,) x Gi"\c). (4.20) 

l^j<fc^2n j=l 

Here is a constant which will be determined below, and the scalar function gb{z) has the properties 

g,{z) = Mz) = (4.21) 

The explicit form of gb{z) is as follows: 

9b\Z) — / 12 2. ™4 ™8^ 7Z87^^27Z4~~l8l ' l*-^^J 



Thanks to ( 2.22 ) and ( 4.21 ) the first three equations ( 4.15| - pT7|) are rephrased in terms of g["''(C) 



and i?(C) as follows 



Gr(-.. ,o+i,or--)-™-= E i?(o/o+i):^':c:Gr(--- ,o,o+ir--)^^^^^-s (4.23) 

^ — ^ J ' j+i 



fc,„(C2n)Gi"\C',C2n)^'^- = G^^^C' , C2"J )^' , (4-24) 

(x2zr^)^"fc^, (Ci)g1"^ (cr\ c)^^ ^~ = Gi"\:E-2Ci,c)^^^ (4.25) 

Let us define the following t/®^"-valued meromorphic function 

q'C\mo^ E ®••■®^^e.„gi"^(«^lc)^ (4.26) 

ei + -- + e2„=0 

where 

2n 

Qt\w\cY = Q("Hz«ic)^9i"^(^ic), gi"^(^«ic) = n n(^^^- - ^«"')- (4.27) 

Since the factor g^"-'(u'|C) is symmetric with respect to Cj's and WaS, respectively, the _R-matrix sym- 



metry ( 4.23 ) is also valid in the boundary case. 



We wish to find integral formulae of the form 



(4.28) 
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Here, the kernel has the form 



aeAj=l 



where 



z) = i) (-) ij {wz) = ^— j-^- J—. 

Note that ipbiw, z) = ipbiw, z'^^) — ijji,{w^^, z). For the function iJ^"-' (uijC) we assume that 

• it is anti-symmetric and holomorphic in the Wa G C\{x~^r^^}, 

• it is symmetric and meromorphic in the Q G C\{0}, 

• it has the four transformation properties 



dt'\w\x-H,x) 

'd''^\x^^Wa^,w\0 



w„ — rx 



IT 



Wa„ - rx ^^^^ X '-Wa„ ~ XWa 



cii -, _Q II 1-1 Q 



TT ^ 

11 ^-1,„-1 _ 



aeA 



it satisfies the following recursion relation 

(w',a;~iz2„-i|C",C2n-i,ea:^~^C2n-i) 



2n-2 



2n-l H ^ 



1 - rZ2n-l 

Here we also fix the constant c' as follows 



n 



^2n-l 



XWa 



, _ a;2»(a:2;a;4)^(^2.^4^^4)2r. 



(4.29) 



(4.30) 



(4.31) 
(4.32) 

(4.33) 
(4.34) 



(4.35) 



(4.36) 



The function '&l^\w\Q is otherwise arbitrary, and the choice oi 'd'^\'w\() corresponds to that of solutions 



to (y^y^. 

From ( 4.26| , 4.32, 4.34, 4.35| ), we find that the integrand may have poles at 



x^+^'^zj^ 1 sC j s; 2n, fc £ Z^o, 



X 



-5-4fe^-l 



'Z-' 1 J 2n, fc G Z^o, 
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We choose the integration contour Ca with respect to Wa (a G A) such that Ca is along a simple 
closed curve oriented anti-clockwise, and encircles the points x^^^^z^^ (1 ^ j ^ 2n, k e Z^q), but not 

We further fix the contour Ca such that 



(1 j 2n, fc G Z^o) nor x^V"^ 



a; "Zj- * lie outside Ca, as in the bulk model case. 



4.4 Main propositions 

We are now in a position to state the proposition in the boundary anti-ferromagnetic XXZ model case: 
Proposition 2 Assume the properties of the function i?^"^ below (4-3C ) and the integration contour Ca 



below \1.3^ . Then the integral formulae (^.13^ \4.2(\ , \i.28j with ^.23^ \12^, ^.28^ , ^.2% \i.2% \4.3(\ ) 
solves the four equations ( 4-l^ - \4~^^ - 



A proof will be given in Appendix A. 

The following example of i^i"'' satisfies all the properties given below ( |I30| ): 

1?(")(w|C) = n^^" J]^ TT Wbe^2{Wa/wb)Q^2{WaWb) 

j=l a£A 



XrWa 



n 

a, be A 
a<b 



XWaWb 



(4.37) 



Taking into account that our integral variables ui^'s correspond to q^Wa in you can easily confirm 



that our integral formulae with the choice of the function t^^"' (4.37) reproduce the ones obtained in the 
trace construction |12|, up to a trivial constant. 



5 Concluding remarks 



In this paper we have constructed the correlation functions of both the bulk and boundary XXZ model in 
the antiferromagnetic regime. This was done by directly solving Smirnov-type master equations; i.e., the 
i?-matrix symmetry, cyclicity conditions (resp. reflection properties) and the normalization conditions 
for the bulk (resp. boundary) case. Precisely speaking, we made the Ansatz that the correlation 
functions of the bulk (resp. boundary) XXZ model are expressed in terms of integral transform of the 
rational functions Q'^J^\w\C,) (resp. Q["''(w|C)) with the integral kernels ^'i"'' (wjC) (resp. '^^^\w\C,)). 
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The integral solutions to Smirnov-type master equations for the bulk case were given by (2.11 
l29| ) with ( ^.23t |2^ , |23o| , |3|), where the function 7?l"^(w|C) satisfies ( |2.32H2.34| ). Those for the 



boundary case were given by ( [4.13t ^4.20| , |^) with ^4.22| , |2.28| , |4.27| , [4.29| , [4.30D , where the 
function (w|C) satisfies ( 1.31 - 1.35 ). Our integral formulae with appropriate choice of the functions 



'!?i"^(w|C) E^nd t?b"''(wlC) reproduce the known results O, |2|. The explicit choice was given in (3.10) 



(resp. ( 1.37 )) for the bulk (resp. boundary) case. 

In comparison with the form factor integral formulae ^, |^] , the corresponding rational function has 
the determinant structure and consequently different from Q^(P'\w\C,) obtained in the present paper. On 
the other hand, the integral kernels are essentially the same. Now we have a natural question whether 
the integral formulae for form factors of the boundary XXZ model can be written by using the integral 
kernel '^^^\w\C,). We wish to address this problem in a subsequent paper. 
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A Appendix A. A proof of Proposition |2| 



The i?-matrix symmetry ( 4.15| ) can be shown by repeating the similar arguments as proving (2.13) 



Straightforward calculations shows that the normalization condition (4.18) follows from (4.35, 4.36) 
Here the formulae (3.4-3.7) are again useful. In order to show the reflection properties ( 4.16| , 4.17) 



the following observation is very relevant. The recursion relation ( 4.35 ) indicates that the function -i^^"-* 
contains the factor 



aeA 



xrwa 



Thus we find that 



i?["''(w|C) has zeros at Wa = x^^ [n G Z), and has a pole at Wa 



(A.l) 



Let us prove ( |4.16D by showing (4.24). First consider the case £2n < 0. By comparing the integrands 



of both sides, the desired equality follows from (4.31) and the equality 



'dt'\MCX2n) WixWa - Z2n){xZ2n - lU^ ^ ) = , C2n) W^ixWa 



^2n ) {^^2n 



When £2,1 > the relation (4.24) is equivalent to 

Z2n ~ r 

The difference between both sides are equal to 



-n 



2'Kiw, 



1 - rZ2n 

1 - rz:;} 



XW2n - Z2n 
XW2n - ZnJ- 



(A.2) 



Since the integral (A.2) vanishes, the relation (4.24) with £2ri > holds. In order to see this, let us 
consider the change of the integral variable W2n '-^ ^271 ■ Under this transform, the integral contour 



C2n is invariant while the integrand of (A.2) negates the sign. The latter follows from (4.31) and the 
equality 



E <"^K'^2jC',C2n)(^«27 - n i^^'^L ^^a)-' = 0. 

Thus the one of the reflection properties ( 4.16|) is proved. 



aeA 



Finally, let us prove another reflection property (4.17) by showing (4.25). First consider the case 



£i < 0. Since the function Qf^ {w\C) contains the factor {xzi — 'Wa){xzi ~ for any a Cz A, there 

is no pole at Wa = (xzi)^^ in the integrand of g["^(C)~^. When the RHS of ( [4.25| ) with £i < is 
analytically continued from (^i to the poles of the integrands 



„13±1 9±1 ^5,±1 



-9±1 



, X Z-y , X Z-^ •) X Z^ , XZ-y ^ X Z\ ^ X Z-y ^ X Z-y , 



-13„±1 



18 



move to the points 



„13^±l ^9^±1 _5,±1 
'1 ' 



X Zj^ ^ X Z-^ ^ X Z-^ ^ X Z\ ^ X Zj^ ^ X Z-^ ^ X Z-^ ^ 

Thus the integral contour Ca for any a £ ^ is invariant even if we make the analytical continuation 



Ci 



^Ci. Furthermore, by using (4.33), the desired equality follows from that of the integrands. 



When £i > the LHS of (|^) times (x^^zi)^"-! reduces to 



n 



dwg 
2TriWn 



r(-icr\c)Qr'(-icr\c)Y 



1 — X rz^ 



(A.3) 



Concerning the RHS of (4.25), the integral contour Ca for a S ^\{1} is again invariant from the same 
reason as the previous case. As for the the integral with respect to wi, we have to keep in mind that 
the integrand of G^jf' {C,)^^ has a pole at wi = xzi. Thus the RHS of (j^H) times (x-^zi)^"-! reduces 
to 

^vi;(")(^icr\c)Qi"^(^lcr^-'^^^ 



dWa 



-27riRes„,j=^-32^ 



Wi 



X •^zi — Wi 



(A.4) 



Let us denote (A^) by L and ( A.4 ) by R, respectively. Furthermore, we will divide R into the sum of 
three parts as follows. In (A.4) we denote the integral part with respect to wi by and the residue 
part by i?2 Rs, the difference of two contour integrals. For that purpose, we introduce the function 



R{w) 



5 —4 —1 1 

zi — X Wi 1 — xrwi 



Wi - x^^w^ ^ I — X "^rzi ' 

as we have done in Q. The integrand of (A.4) has poles at wi = x~^zi,xzi^ ,x^^r^^ in the annulus 
between the closed lines Ci and x^'^Ci. Here, the closed line a;"'*Ci denotes the one obtained from Ci by 
similarity transformation with the ratio being x~^. Since R{x-^zi) = 1 and i?(a;zf ^) = R{x-^r-^) = 0, 
the residue at wi = x^^zi can be evaluated as follows: 



n 

aeA\{l} 



dWa 



dwi ^^^(n), 

- b 



X ■^Zi — Wi 



R{wi). (A.5) 



Let us denote the first and the second term of ( [A.5| ) by i?2 and — i?3, respectively. The function R{wi) 
has a pole at wi — x^^ , however, we notice that the integrand of ( |A.5| ) is regular at this point because 



of (A.l). The integral i?2 can be rewritten in terms of the integral along the contour Ci by changing 
the integral variable wi >— > x^^w^^ as follows: 



dwa 



i? = IT 



dwi („) 



_i ~.xz^'^~x '^U'l^ 1 — xrwi 



Thus we obtain 

L — i?i = i?2 ^ R3, 



n 



dWr, 



Wi — x^'^Wi ^ I ~ X ^rzi 
{n)f 1^-1 x^^zi - xzf ^ 1 - xrwi 



X ■^zi — wi 1 — x "^rzi ' 



which implies ( 4.25 ) with ei > 0. Thus Proposition g was proved. 
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